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DEBRIS  DISTRIBUTION 
SUMMARY 

INTRODUCTION 

This  study  is  a  continuation  of  the  ’’DEBRIS  CLEARANCE  STUDY.  "* 
That  study  developed  methods  of  estimating  gross  debris  accumulation  in 
various  types  of  urban  areas,  based  on  total  structural  demolition  and  uni¬ 
form  distribution  of  the  debris  over  a  given  size  area. 

The  present  study  intends  to  improve  the  degree  of  sophistication  of 
predicting  debris  distribution.  The  modus  operandi  of  this  improvement 
was  to  divide  debris  into  two  categories;  stationary  and  transportable. 
Application  of  fragmentation  and  trajectory  models  was  made  on  the  trans¬ 
portable  debris  and  distribution  was  based  on  this  type  of  debris.  Although 
transportable  debris  may  originate  from  such  sources  as  buildings,  auto¬ 
mobiles,  communication  equipment,  public  utilities,  and  any  other  exposed 
object,  this  report  addresses  itself  exclusively  to  buildings.  To  be  specific, 
brittle  wall  elements  of  buildings  were  studied  in  detail.  Because  of  blast 
orientation,  thermal  effects  and  material  properties,  brittle  wall  elements 
were  considered  to  be  one  of  the  primary  sources  of  transportable  debris. 

A  descriptive  outline  for  estimating  debris  distribution,  in  an  area 
of  interest,  was  developed.  It  was  shown  that  if  one  can  specify  the  relation¬ 
ship  between  load  and  fragment  piece  size  distribution  (i.  e.  ,  fragmentation) 
then  area-wide  debris  distributions  may  be  developed  with  a  fair  degree  of 
confidence. 

In  order  to  predict  the  fragmentation  of  a  brittle  wall  element  an 
analytical  model  was  developed.  This  model  was  based  on  a  probabilistic 
approach  and  the  basic  underlying  assumptions  were  verified  by  experi¬ 
mentation.  Debris  profiles  of  transportable  debris  were  developed  utili¬ 
zing  *.he  analytical  fragmentation  model  and  a  trajectory  analysis. 

*Ahlers,  E.  B.  ,  "Debris  Clearance  Study,  ”  IIT  Research  Institute,  for  the 
Office  of  Civil  Defense,  OCD-FS-62-202,  Subtask  3322-A,  September  1962. 
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In  order  to  demonstrate  the  feasibility  of  the  techniques  developed 
a  hypothetical  structure  was  analyzed.  This  structure  consisted  of  a  single 
wall  composed  of  similar  masonry  wall  panels.  The  individual  panels  were 
8  ft- 6  in.  high  and  6  in.  thick.  This  panel  when  suitably  idealized  and 
uniformly  loaded  yielded  the  fragmentation  characteristics  as  shown  in 
Fig.  S.  I.  This  curve  illustrates  that  fragmentation  occurs  at  about  8  poi 
and  total  fragmentation  is  achieved  at  about  20  psi.  Figure  S.  2  illustrates 
the  transport  characteristics  of  the  resulting  panel  fragments  under  similar 
loading  conditions  and  at  various  heights  above  ground.  Information  obtained 
from  Fig.  S.  1  and  S.  2,  together  with  some  knowledge  of  wall  dimensions 
are  sufficient  to  develop  a  debris  profile  for  the  entire  wall  as  shown  in 
Fig.  S.  3.  Walls  of  various  heights  and  fragmentation  characteristics  were 
investigated  and  the  results  summarized  in  Table  S.  1. 

Table  S.  1 

SENSITIVITY  OF  DEBRIS  CONTOURS  TO  BUILDING  HEIGHT  AND 

PARTICLE  SIZE 


Building 

Height 

All  Small 

Predicted 

All 

Large 

* 

X 

.  y 

# 

X 

* 

y 

* 

X 

* 

y 

40 

514 

0.  003 

200 

0.  036 

200 

0.  085 

35 

496 

0.  0027 

180 

0.  035 

180 

0.  080 

30 

460 

0. 0023 

1  65 

0.  032 

130 

0.  075 

25 

425 

0.  019 

155 

0.  029 

— 

— 

20 

400 

0.  015 

135 

0.  026 

130 

0.  060 

15 

350 

0.  012 

125 

0.  023 

125 

0.  055 

10 

300 

0.  009 

75 

0.  018 

90 

0.  045 

5 

175 

0.  006 

50 

0.  015 

40 

0.  040 

x  .  y  are  the  coordinates  of  the  peak  value  of  the  debris  profile. 
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FRAGMENTATION  CURVE  FOR  MASONRY  PANEL 


Fig.  S.  2  FRAGMENTATION  AND  TRANSPORT  CHAR ACTERISTICS  FOR  MODEL  PANEL 


The  findings  of  this  study  may  be  summarized  as  follows. 

1.  Debris  may  be  divided  into  two  categories,  offsite  and  onsite. 
Onsite  debris  prevents  no  impedance  to  logistic  activity  other  than  localized 
shelter  rescue  operations. 

2.  In  order  to  characterize  offsite  debris  profiles,  the  piecesize 
distribution  of  the  fragments  must  be  known.  This  is  in  c  rder  to  character¬ 
ize  the  particle's  trajectory. 

3.  Low  level  dwellings  must  be  included  as  potential  debris  pro¬ 
ducing  structures.  Even  elements  at  10  ft  above  ground  may  be  transported 
significant  distances  (e.  g.  ,  into  the  street)  at  moderate  overpressure 
levels  (10  psi). 

4.  In  general,  smaller  particles  from  higher  initial  heights  will 
be  transported  greater  distances  than  larger  particles  at  lower  heights. 

As  overpressure  increases  the  size  of  the  particle  fragmenting  will  be 
smaller. 

5.  Particles  fragmenting  and  being  transported  offsite  from  their 
original  position  will  have  an  extremely  high  terminal  velocity  (i.  e.  ,  at 
least  50  ft  per  sec).  For  this  reason,  offsite  debris  must  be  considered 
an  extremely  dangerous  secondary  effect  of  blast.  This  effect  must  be 
evaluated  in  light  of  its  damage  producing  capabilities  to  structures  as 
well  as  a  casualty  producing  mechanism  to  unsheltered  populenca. 

6.  The  assumption  that  the  maximum  dynamic  stresses  introduced 
into  the  various  unit  segments  are  independent  of  the  fracture  characteristics 
has  been  shown  both  experimentally  and  analytically.  This  establishes 

the  theoretical  development  of  the  fragmentation  model. 

7.  The  statistical  strength  characteristics  of  typical  unit  segments 
have  been  expressed  by  their  cumulative  distribution  function,  F(p).  This 
function  gives  the  probability  of  fracturing  the  unit  at  a  load  magnitude 
equal  to  or  less  than  p  and  for  computational  purposes  has  been  expressed 
in  Weibull  form.  It  is  tc  be  emphasised  that  this  form  was  chosen  for  con¬ 
venience.  At  present,  Weibull  parameters  are  not  Known  for  many  common 
construction  materials  and  must  be  determined. 
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8.  Debris  profiles  from  multistory  walls  may  be  generated.  The 
profiles  are  typically  unimodal  in  shape  and  of  course  are  very  sensitive 
to  piecesize  distribution.  For  walls  at  5,  10,  15  stories  in  height  profiles 
exist  for  each  of  three  piecesize  distributions,  all  large,  predicted,  and 
all  small.  It  is  interesting  to  note  that  the  large  particles  generally  lend 
most  of  the  shape  to  the  predicted  size  profile.  Thai  is,  up  to  the  peak  the 
large  size  particles  predominate.  This,  of  course,  is  in  part  due  to  the 
method  of  developing  these  profiles.  Large  particles  are  distributed  over 
much  less  area  than  small  ones  and  hence  tend  to  produce  much  greater 
depths. 

9.  The  profiles  have  their  ordinate  normalized  with  respect  to 
volume  and  packing.  (Packing  corresponds  to  a  void  ratio  of  l.  0).  If  the 
ordinate  were  multiplied  by  the  unit  width  volume  of  material  at  one  story 
height  and  then  again  by  a  suitable  void  ratio,  then  the  profile  would  ex¬ 
press  the  true  depth  of  the  building's  transported  debris  at  all  points  along 
its  transported  distance. 

The  analytical  methods  of  fragmentation  and  subsequent  transport 
did  not  consider  either  the  effects  of  orientation  of  the  structure  to  the 
blast  or  shielding  of  one  building  by  another.  It  is  therefore  recommended 
that  any  follow  on  study  should: 

1.  Continue  to  develop  analytical  methods  for  additional 
structural  materials  and  elements  as  well  as  complete 
structures.  These  analytical  models  should  be  suitably 
verified  by  appropriate  experimental  investigation. 

2.  The  individual  building  deoris  contours  mast  be  com¬ 
bined  to  give  a  cumulative  debrie  contour  for  an 
entire  subarea  of  contiguous  structures. 

3.  Effects  of  blast  orientation  and  structural  shielding 
should  be  accounted  for  in  specifying  the  load  on  the 
structure. 

4.  Finally,  all  methods  developed  ehould  be  combined  to 
yield  a  single  computational  model. 

lit  aiSSAtCM  INI! nuts 

S-7 


DEBRIS  DISTRIBUTION 


ABSTRACT 

This  study  is  a  continuation  of  the  "DEBRIS  CLEARANCE  STUDY."1* 
That  study  developed  methods  of  estimating  gross  debris  accumulation  in 
various  types  of  urban  areas,  based  on  total  structural  demolition  and 
uniform  distribution  of  the  debris  over  a  given  size  area. 

The  present  study  makes  use  of  several  of  the  tools  developed  in 
the  previous  report  and  expanded  herein  to  more  accurately  predict  the 
debris  distribution.  Structural  fragmentation  is  defined  for  a  particular 
set  of  blast  and  structural  parameters.  The  trajectories  of  the  transport¬ 
able  material  are  calculated,  and  the  distribution  of  this  material  is  found 
by  size  as  well  as  mass.  In  order  to  establish  the  credibility  of  the  frag¬ 
mentation  theories  which  were  developed,  a  series  of  experiments  were 
conducted  on  brittle  beams  under  dynamic  loading. 

A  sample  problem  is  worked  out  that  utilizes  the  above  techniques 
on  a  hypothetical  building  to  illustrate  how  the  building's  debris  distribution 
may  appear  under  various  conditions  of  building  height  and  fragment  size 
distribution. 
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INTRODUCTION  AND  ORGANIZATION 


If  there  is  to  be  any  simulation  of  the  postattack  environment  of  any 
given  locale,  certain  input  information  must  be  available  to  the  simulation 
model.  Among  this  input  may  be  such  things  as  specifying  an  urban  area 
analytically,  defining  clearance  priorities,  and  specifying  debris 
distributions. 

This  report,  although  it  is  a  continuation  of  the  previous  postattack 
report:  "DEBRIS  CLEARANCE  STUDY,  intends  to  improve  the  degree  of 
sophistication  of  predicting  debris  distribution.  The  modus  operand!  of 
this  improvement  was  to  divide  debris  into  two  categories:  stationary  and 
transportable.  Application  of  the  fragmentation  and  trajectory  models  was 
made  on  the  transportable  debris  and  distribution  will  be  based  on  this 
type  of  debris. 

As  an  example  of  the  application  of  the  techniques  developed,  a 
hypothetical  exterior  masonry  wall  was  represented  by  an  equivalent  beam 
analogy  and  debris  profiles  computed  for  various  heights  and  size  distribu¬ 
tions. 

The  report  is  organized  into  five  chapters.  Chapter  One  is  a  des¬ 
criptive  outline  or  plan  of  attack  for  estimating  debris  distribution  in  a  given 
aiea  of  interest.  It  will  be  shown  here  that  if  one  can  specify  the  relation¬ 
ship  between  load  and  fragment  piecesize  distribution  (i.  e.  ,  fragmentation; 
then  area-wide  debris  distributions  may  be  developed  with  a  fair  degree 
of  confidence. 

The  phenomena  of  structural  fragmentation  of  brittle  materials 
(e.  g.  ,  masonry  and  concrete)  is  discussed  in  Chapter  Two.  The  fragmenta¬ 
tion  model  for  various  brittle  beam  conditions  is  presented  and  the  applica¬ 
tion  of  its  results  are  discussed.  Review  of  the  transport  model  is  covered 
here  too. 

Chapter  Three  attempts  to  establish  the  validity  of  the  fragmentation 
theory  dibcussed  above.  Brittle  beams  with  various  end  conditions  were 
subjected  to  a  variety  of  external  dynamic  loadings. 
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An  example  of  a  hypothetical  masonry  wall  is  investigated,  utilizing 
the  methods  of  the  previous  chapters,  in  Chapter  Four.  Also  developed  is 
a  procedure  which  allows  one  to  develop  debris  profiles. 

Finally,  Chapter  Five  contains  a  summary,  conclusions  and  recom¬ 
mendations  concerning  the  work  done  in  this  report. 
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CHAPTER  ONE 


A  MODEL  FOR  PREDICTING  DEBRIS  DISTRIBUTION 
1.  GENERAL 

In  order  to  predict  the  distribution  of  offsite  debris  in  a  metropolitan 
area,  a  great  deal  of  information  must  be  known.  Lack  of  this  information 
does  not  preclude  the  possibility  of  using  intuitive  assumptions.  These 
assumptions,  of  course,  must  be  pointed  out  as  such  and  some  qualitative  , 
if  not  quantitative,  check  should  be  made  on  them. 

In  this  report  we  wish  to  investigate  the  quantity  and  composition 
of  offsite  debris  (e.  g.  ,  debris  which  is  transported  from  its  original  loca¬ 
tion)  which  hampers  postattack  logistical  missions.  Although  such  debris 
may  originate  from  such  sources  as  buildings,  automobiles,  communication 
equipment,  public  utilities,  and  any  other  exposed  object,  this  report  will 
address  itself  exclusively  to  buildings.  To  be  specific,  brittle  wall  elements 
of  buildings  will  be  studied  in  detail.  Because  of  blast  orientation,  thermal 
effects,  and  material  properties,  brittle  wall  elements  are  considered  to 
be  one  of  the  primary  sources  of  the  transportable  debris. 

The  following  sections  of  this  chapter  describe  a  step-by-step  ap¬ 
proach  by  which  a  postattack  analyst  might  map  the  debris  distributions  for 
an  entire  metropolitan  area.  This  approach  is  not  to  be  considered  unique. 
Any  other  rational  method  of  analysis  is  not  precluded. 

L.  1  Division  of  Metropolitan  Areas  by  Representative  Building  Types 

By  use  of  aerial  photographs,  Sandborne  maps,  or  some  other  con¬ 
venient  source  the  investigator  will  find  that  most  metropolitan  areas  may 
be  broken  down  into  sub-areas  which  may  be  represented  by  one  or,  at 
most,  only  a  few  different  types  of  structures.  Such  a  sub-area  may  be 
considered  homogeneous  from  the  standpoint  of  potential  debris  production. 

It  now  becomes  necessary  to  investigate  the  different  types  of  building 
construction  in  order  to  evaluate  the  condition  of  each  sub-area. 

1.  I  Specification  of  Blast  Loading  on  a  Single  Building  Type 

When  one  investigates  the  formulation  and  transportation  of  debris 
from  a  single  source  the  loading,  due  to  primary  weapon  effects,  must  be 
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available  in  a  manner  such  as  the  following: 

P  =  P„  (W,  X  h.;  k) 

O  O  O  D 

where, 

Pq  =  Peak  Overpressure 

W  =  Weapon  Yield 

X  =  Distance  from  Ground  Zero 
o 

h^  =  Height  of  Burst 
k  =  Shielding  Factor 

The  description  of  the  above  relationship  has  been  the  subject  of  many 
2  3 

investigators.  ’  One  of  the  less  highly  investigated  areas  is  the  relation 

ship  between  the  shielding  factor,  k,  and  the  degree  of  structural  density 

at  the  distance  X  . 

o 

1.  3  Debris  Production  for  a  Single  Building 

The  next  step  in  establishing  debris  levels  for  structures  is  to 
describe  the  actual  distribution  of  material  in  a  particular  building.  This 
may  be  shown  in  the  manner  illustrated  in  Fig.  1.  1.  Each  curve  is  sym¬ 
bolic  of  different  types  of  materials  and  might  include  concrete,  masonry, 
plaster,  etc. 

It  becomes  necessary,  next,  to  bo  able  to  predict  the  piecesize 
distribution  of  these  materials  as  a  function  of  the  loading  on  the  structure. 
Such  a  relationship  is  shown  in  Fig.  1.  2. 

1.  1  Final  Deposition  of  Debris  Material  for  a  Single  Building 

If  the  initial  position  and  size  of  debris  is  known,  it  is  possible  to 
accurately  predict  its  transported  position  on  the  ground  by  computing  its 
trajectory.  In  order  to  get  a  picture  of  how  the  debris  accumulates  or 
distributes  something  must  be  known  about  the  available  open  area  sur¬ 
rounding  the  debris  particle's  final  transported  position.  That  is,  in  an 
open  region  debris  can  spread  out  and  not  be  as  deep  as  it  might  be  under 
more  crowded  conditions. 

It  is  possible,  therefore,  to  draw  the  distribution  curve  of  offsite 
debris,  for  a  single  structure,  if  a  one  to  one  correspondence  between 
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debris  produced  and  final  ground  position  is  known.  Such  a  curve  may 
very  well  turn  out  to  look  like  the  one  shown  in  Fig.  1.  3.  It  is  observed 
that  the  maximum  debris  level,  y  ,  occurs  at  some  distance,  x  ,  from 
the  structure  for  a  particular  loading. 

1.  5  Estimating  Level  of  Constant  Debris  for  Sub-Area 

Now  that  debris  profiles,  as  shown  in  Fig.  1.  3,  are  available  for 
each  building  type  in  a  particular  sub-area,  a  single  cumulative  debris 
height  may  be  determined  from  the  envelope  of  the  superimposed  individual 
profiles.  Such  a  relationship  is  showm  in  Fig.  1.4. 

If  x _  represents  the  maximum  debris  level  distance  for  an  indi¬ 

max 

vidual  debris  profile  making  up  the  cumulative  curve,  then  this  distance 
will  represent  the  inner  boundary  of  the  original  sub-area.  Within  the  inner 
sub-area  there  will  be  a  constant  level  of  debris  equal  to  the  cumulative 
debris  height,  y  .  This  is  shown  in  Fig.  1.  5  and  is  a  conservative 
estimate.  Figure  1.  6  represents  these  "islands"  of  constant  debris  height 
within  the  overall  area.  It  remains  to  interpolate  the  debris  depths  between 
the  islands. 

1.  6  Interpolation  Based  on  Blast  Angle  of  Incidence 

Up  to  this  point  it  has  been  assumed  that  the  blast  impinges  at 
normal  incidence  on  the  structure.  This  assumption  was  made  to  simplify 
the  loading  analysis.  At  this  point,  however,  the  blast  angle  of  incidence 
(i.  c.  ,  to  each  "debris  island")  may  now  be  introduced  to  help  weight  the 
process  of  interpolation  of  the  debris  level  between  islands.  The  final 
result  of  such  a  process  leads  to  a  map  of  debris  level  for  the  entire  area 
under  investigation.  Such  a  map  is  pictured  in  a  qualitative  manner  in 
Fig.  1.7. 
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cumulative:  offstte  debris  profile  for  a  sub-area 


.  S  FORMATION  OF  "DEBRIS  ISLAND"  WITHIN  A  SUB-AREA 


CHAPTER  TWO 


AN  ANALYTICAL  APPROACH  TO  PREDICTING  SIZE 
AND  TRANSPORT  OF  OFFSITE  DEBRIS 

2.  GENERAL  APPROACH 

This  chapter  deals  with  a  problem  that  is  fundamental  in  any  analyti¬ 
cal  model  of  i .ie  postattack  environment.  As  has  been  discussed  in  the  first 
chapter  of  this  report,  one  of  the  most  predominant  forms  of  offsite  debris 
will  be  masonry  and  other  forms  of  frangible  rubble.  In  order  to  predict  the 
final  deposition  of  known  sources  of  this  offsite  debris  the  actual  mode  of 
fracture  must  be  investigated  as  well  as  the  subsequent  motion  of  the  frag¬ 
mented  element. 

As  has  been  outlined  previously*  the  prediction  of  the  formation 
and  transport  of  offsite  debris  is  a  four  step  process  that  includes: 

1)  The  prediction  of  the  reflected  and  drag  pressures 
resulting  from  a  given  set  of  attack  parameters. 

2}  Determination  of  the  reflected  and  drag  pressures 
on  various  structural  elements. 

3)  The  prediction  of  when  and  how  a  brittle  structural 
element  will  fail. 

4)  The  motion  of  the  displaced  element  in  space  unrit*r 
gravitational  and  wind  forces. 

The  problem  posed  by  the  first  two  steps  of  the  description  above 
has  been  taken  up  satisfactorily  in  a  number  of  well  known  documents.^’  ’  * 
The  final  two  steps  are  of  a  more  fundamental  nature  anti  will  be  dealt  with 
here  in  detail. 

Having  modeled  the  actual  formation  and  transport  of  offsite  debris 
it  becomes  evident  that  a  one-to-one  correspondence  can  be  established 
between  materials  in  tin  target  and  their  final  position  on  the  ground  surface. 
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2.  1  A  Simple  Fragmentation  Model 
2.  1.  1  Introduction 


The  ultimate  strengths  obtained  from  the  repeated  static  testing  of 
nominally  identical  brittle  specimens  will  exhibit  a  characteristic  scatter. 
Furthermore,  the  locations  of  the  resulting  fractures  will  vary  from  one 
specimen  to  another.  Because  the  disposition  of  ultimate  strengths  is 
usually  very  large  for  brittle  materials,  it  is  generally  not  possible  to 
predict  the  behavior  of  a  single  element  with  any  useful  accuracy  regardless 
of  the  amount  of  accumulated  experience  with  similar  elements.  It  is  pos¬ 
sible,  however,  through  the  use  of  statistics,  to  predict  the  composite 
behavior  of  a  large  group  of  nominally  identical  brittle  components  from  a 
knowledge  of  the  characteristics  of  still  another  large  group  of  similar 
components. 


Adopting  a  statistical  viewpoint,  we  shall  consider  the  fragmentation 
of  a  brittle  beam  structure  such  as  the  cantilever  shown  in  Fig.  Z.  1.  Under 
a  sufficiently  high  static  loading,  the  statistical  theory  of  fracture  allows 
that  fracture  may  occur  anywhere  in  the  span.  The  likelihood  of  fracture 
will  be  greatest  at  the  fixed  end  where  the  stress  level  is  highest.  It  can 
be  observed  that  the  probability  of  fracture  <  ccurring  at  a  specified  station 
along  the  beam  is  zero.  Any  finite  probability  of  failure  at  the  various  beam 
stations  would  always  result  in  the  physical  contradiction  that  the  survival 
probability  of  the  beam  is  zero.  Consequently,  it  is  meaningless  to  seek 
fragments  of  a  given  size;  and  indeed,  the  number  of  fragments  which  occur 
within  a  given  range  of  sizes  should  be  sought. 
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2.  1.2  Method  of  Attack 


To  obtain  a  finite  probability  of  fracture,  we  must  consider  a  region 

along  the  beam  rather  than  a  specific  station.  Therefore,  we  shall  imagine 

the  beam  to  be  divided  into  n  units  of  equal  length.  T<  each  of  these  units 

will  be  assigned  a  measure  of  its  statistical  strength  characteristics.  This 

will  take  the  form  of  a  cumulative  distribution  function,  F  (p),  which  is  a 

measure  of  the  probability  of  fracturing  the  unit  at  a  load  magnitude  equal 

or  less  than  a  load  p.  It  remains  to  investigate  the  probability  of  obtaining 

any  size  element  resulting  from  single  or  multiple  fractures  occurring  within 

any  of  the  n  subunits  of  the  beam.  It  will  be  seen  that  the  problem  of  pre- 

5 

dieting  fracture  is  an  analogy  of  the  classical  "run"  probability  problem. 

2.1.3  Dynamic  Stresses 

The  fundamental  assumption  in  the  development  of  the  fragmentation 
model  is  that  the  maximum  dynamic  stresses  introduced  into  the  various 
unit  segments  are  independent  of  the  fracture  characteristics  of  these 
segments.  This  would,  of  course,  be  the  case  if  the  loading  was  sufficiently 
rapid  to  fully  stress  the  unit  before  any  fracture  occurred.  This  important 
assumption  leads  to  the  simplification  that  the  units  can  be  treated  as  being 
stochastically  independent.  For  static  loading,  this  assumption  is  clearly 
invalid  since  fracture  in  a  unit  segment  would  immediately  relieve  the 
stresses  in  other  segments.  This  implies  that  only  one  fracture  can  occur 
in  a  statically  loaded  statically  determinate  beam.  On  the  other  hand,  even 
the  crudest  dynamic  loading  experiments  with  statically  determinate  brittle 
beams  produce  multiple  fracture. 

Taking  the  dynamic  beam  h  ading  in  the  form' 

P  f  (x)  g  (t) 

where  the  effects  of  magnitud  .  loud  distribution,  and  pulse  shape  are  ex- 
plit  itly  delineated  the  maximal  i  dynamic  stresses  c.m  be  estimated  lor  a 
typical  unit  !>'  multiplying  the  static  stresses  associated  with  p  f  (x)  by  the 
maximum  dynamic  load  factor  (DLF)  for  the  first  mode  >f  vibration.  Thus, 
tin-  stress  distribution  f<  r  a  given  unit  segment  ha>  the  form 

n  h  (x.  v.  /)  (DLF)  (See  Appendix  A) 
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2.  1.4  Statistical  Strength  of  Unit  Segments 

The  statistical  strength  characteristics  of  a  typical  unit  segment 
are  described  by  its  cumulative  distribution  function,  F  (p).  This  function 
gives  the  probability  of  fracturing  the  unit  at  a  load  magnitude  equal  to  or 
less  than  p.  The  distribution  function  may  be  determined  by  physically  testing 
many  unit  segments  or  by  appealing  to  one  of  the  "weakest  link"  statistical 
fracture  theories.  Of  these,  Weibull's6  is  the  most  famous  theory  and  is 
shown  below  in  a  form  which  is  appropriate  for  the  unit  segment. 


o’-  <r 


m 


u 


d  V 


F  (p)  =  1  -  e' 
=  0 


o 


cr  >  a 


(2-1) 


u 

ru 


where  is  the  stress  distribution  in  the  segment  of  the  form  ph  (x,  y,  z)(DLF); 

<r  <r  and  m  are  constants  of  the  material;  and  where  the  integration  is 
u  o 

taken  over  the  volume  of  the  unit  segment.  The  distribution  parameters 

<r  ,  o-  ,  and  m  are  usually  determined  from  simple  bending  or  tension 
u  o 

specimens.  (Appendix  B) 


The  practical  application  of  the  Weibull  function  takes  advantage  of 
certain  mechanical  properties  normally  found  in  most  brittle  materials. 

For  example,  the  linearity  of  the  stress  -  strain  relationships  up  to  rupture 
greatly  simplifies  the  stress  analysis.  Also  the  insensitivity  of  the  strength 
of  the  brittle  materials  to  wide  variations  in  strain-rate  enables  the  deter¬ 
mination  of  a  unique  set  of  distribution  parameters  from  static  test  results. 
Finally,  since  the  tensile  strength  is  generally  much  smaller  than  the  com¬ 
pressive  or  shear  strength,  a  consideration  of  tensile  bending  stresses  alone 
is  usually  sufficient  for  the  determination  of  the  strength  of  normally  pro¬ 
portioned  beams. 


/ 


2.  1.  5  The  Problem  of  Runs 

When  a  beam  undergoes  multiple  fracture  along  its  span,  different 
size  elements  occur  due  to  the  proper  combination  of  fractures  taking  place. 
If  one  were  to  examine  each  combination  of  possible  fracture  this  would 
involve  2°  -  1  possibilities.  The  result  of  the  investigation  would  lead  to 
being  able  to  predict  the  probability  of  getting  an  element,  of  a  given  size 
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range,  from  a  particular  part  of  the  loaded  beam.  If  it  is  sufficient  to 
only  know  the  probability  of  getting  an  element,  with  a  given  size  range, 
from  an  unspecified  point  along  the  original  beam,  the  fragmentation  problem 
is  directly  analogous  to  the  problem  of  runs. 

A  run  of  exactly  r  successes  will  be  said  to  occur  if  it  is  possible  to 
find  a  sequence  of  r  consecutive  successes  which  is  not  part  of  a  sequence 
■of  more  than  r  successes.  Thus  a  run  of  exactly  r  successes  can  occur  either 
at  the  beginning  of  the  scries  in  which  case  it  must  be  followed  by  a  failure, 
or  at  the  end  in  which  case  it  must  be  preceded  by  a  failure,  or  elsewhere 
in  which  case  it  must  be  both  preceded  and  followed  by  a  failure. 

No  difficulty  arises  in  calculating  the  probability  that  a  run  of  exactly 
r  successes  will  occur  in  n  trials  so  long  as  r^l/Z  n,  (i.e.,  so  long  as  it 
is  not  possible  for  two  or  more  separate  runs  to  occur  in  the  same  series). 
For  example,  a  run  of  exactly  five  successes  in  eight  trials  can  occur  as 
follows,  where  p  denotes  success,  q  denotes  failure  and  X  denotes  an  un¬ 
specified  result: 

5 

p  p  p  p  p  q  XX,  probability  p  q  , 

5 

qpppppq  X,  probability  q  p  q  , 

5 

X  qpppppq,  probability  q  p  q  , 

5 

XXqppppp,  probability  q  p 

These  four  ways  are  mutually  exclusive  and  exhaustive  and  the  required 

5 

probability  is  thus  the  sum  of  the  four  separate  probabilities  (i.  e.  ,  2  p  q 
(  1  +  q)).  In  general,  if  r 1/2  n,  the  probability  of  a  run  of  exactly  r  suc¬ 
cesses  beginning  at  the  initiation  of  the  series  or  ending  at  the  end  of  the 

series  is  p  q  in  each  case,  while  the  probability  of  such  a  run  elsewhere 

r  Z 

in  the  series  in  a  specified  position  is  p  q  .  As  there  are  (n  -  r  -  1)  pos¬ 
sible  positions  apart  from  the  beginning  and  end,  the  probability  p  (r)  say 
of  a  run  of  exactly  r  successes  is  given  by 

p  (r)  =  2  pr  q  +  (n  -  r  -  1)  pf  q"  (2-2) 

provided  1/2  n^r^n.  If  r  =  n,  p  (n)  =  pn  obviously. 

It  is  easy  to  see  that  Eq.  (2-2)  always  gives  the  expected  number  of 
runs  of  exactly  r  successes  in  a  scries  of  n  trials,  whether  or  not  r^  1/2  n. 
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First,  suppose  that  1/2  n;  then  the  probability  that  there  is  exactly  1  run 
is  given  by  Eq.  (2-2),  while  the  probability  of  any  higher  number  of  runs  is 
zero,  so  that  Eq.  (2-2)  is  the  expected  number  of  runs.  Suppose  that  r<C  1/2  n 
so  that  more  than  one  run  of  exactly  r  successes  may  occur  in  a  series  of  n 
trials.  Then  in  a  series,  the  number  of  runs  beginning  at  the  first  trial  is 
obviously  either  0  or  1,  and  the  probability  that  it  is  1  is  p  q.  Similarly,  the 

Y 

expected  number  of  runs  beginning  at  the  first  trial  is  p  q.  Similarly,  the 
expected  number  beginning  at  the  second  trial  is  q  p  q,  as  it  is  for  the  third, 
fourth,  ....  (n  -  r  +  1)  th  trial  (i.  e.  ,  a  run  at  the  end  of  the  series  is  p  q). 
Thus  the  expected  number  of  runs  in  the  series  is  Eq.  (2-2).  Of  course, 
some  of  these  possibilities  are  mutually  exclusive,  but  this  is  immaterial 
since  we  are  finding  the  expected  number  of  runs  in  the  series  and  not 
the  probability  that  a  run  will  occur. 

For  the  purposes  of  analysis,  it  is  assumed  that  the  beam  is  sub¬ 
divided  into  n  equal  segments  and  that  cracks  can  occur  only  at  the  n  +  1 
joints  shown  in  Fig.  2.  2. 


i  -s 

n _ i 

n  n  +  1 


Fig.  2.  2  BEAM  NUMBERING  SYSTEM 


The  expected  number  of  fragments,  r  segments  in  length,  is  equal  to  the 
sum  of  the  expectations  of  such  fragments  starting  at  each  of  the  joints 
1  .  2,  ...  n  -  r  +  1.  That  is, 


whe 


E(r)  = 


n  -  r  +  1 
V 


(2-3) 


Jk=qkVn  f°r  r=1>and 


Jk=  \ 


k  +  r  -  1 


i  I 

i 


) 


[  j=k+i  Ji 


P.  ;  q,  for  r  >  1 

j  i  k  +  r 


(2-4) 


in  which  p.  and  q.  are  the  probabilities  of  the  ith  joint  not  having  or  having 

a  crack,  respectively.  Of  course,  n  ,  .  -  0  and  q  ,  =  1. 

n  +  1  hi  +  1 
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In  the  case  of  a  simple  beam,  Eq.  (2-3)  and  (2-4)  still  hold,  but 

of  course  here  we  have  q  -  1  as  well  as  q  ,  ,=  1.  Where  the  analysis 

is  to  be  made  on  a  fixed-fixed  beam  formulas  (2-3)  and  (2-4)  hold  with  no 

restrictions  on  p  ,  ,  and  q  ,  ,  or  q,. 

rn  +  1  m  +  1  Mr 

The  recursive  relation  (Eq.  2-3)  predicts  the  expected  number  of 
segments  of  each  elemental  size  where  an  elemental  size  ,  r,  consists  of 

r  unit  segments,  (i.  e.  r  =  1,  2,  3 . n  segments  in  length)  Because 

fracture  may  occur  anywhere  inaunit  segment  rather  than  at  the  assumed 
joints  as  shown  in  Fig.  2.  2,  any  elemental  fragment  size  (e.  g.  r  =  1,  2,  3, 

. .  n)  actually  represents  the  midpoint  of  a  range  of  sizes  which  may 

differ  from  the  indicated  size  by  as  much  as  one  unit  length  (  +  ).  It 

is  apparent  that  the  elemental  fragment  size  ranges  overlap  considerably 
as  shown  in  Fig.  2.  3,  and  as  a  consequence,  they  cannot  be  used  directly 
for  classifying  sizes.  To  circumvent  this  problem  we  will  group  the 
elemental  sizes  into  broader  physical  size  intervals  Now',  since  E(r) 
represents  the  expected  number  of  segments  of  elemental  size,  r 
E(r)  x  r  x  L/n  represents  the  total  physical  length  devoted  to  elements  of 
elemental  size  r.  To  convert  this  to  a  fraction  of  the  original  beam  length, 
F(r),  we  divide  by  L,  thus, 

F(r)  =  -£-  E  (r)  (2-5) 

The  expected  fractions  of  the  beam,  as  determined  by  Eq.  (2-5)  may  be 
grouped  into  broader  subdivisions  in  such  a  w'ay  that  the  effect  of  over¬ 
lapping  may  be  made  as  small  as  desired. 

It  is  clear  from  Fig.  2.  3  that  the  relative  amount  of  overlap 
decreases  with  increasing  n.  The  effect  of  overlap  on  the  accuracy  of  the 
fragmentation  predictions  can  be  seen  from  Table  2-1  which  indicates 
little  improvement  after  n  passes  30  trials. 
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Fig.  2.  3  ELEMENTAL  AND  PHYSICAL  SIZE  INTERVALS 


Table  2.  1 


THE  EFFECT  OF  OVERLAP  ON  THE  ACCURACY  OF 
FRAGMENTATION  PREDICTIONS  (P  =  20  psi) 


Size  Range, 
ft 

n  -  5 

15 

30 

50 

100 

6-2 

25.  8 

31.  63 

32.  05 

32.  13 

32.  16 

2-4 

0 

0.  10 

0.  16 

0.  19 

0.  23 

4-6 

14.  6 

6.  87 

4.  58 

3.  77 

3.  22 

6-8 

49.  6 

61.  49 

63.  36 

64.  07 

64.  59 

8-10 

0 

0 

0 

0 

0 

2.  1.6  Hydrostone  Cantilever  Results 

As  a  test  case  of  the  analytical  procedure  described  abo\ 
same  hydrostone  beam  studied  in  the  overlap  considerations  of  the  previous 
section  was  analyzed  for  a  constant  n  (i.  e.  n  =  15)  over  a  complete  range 
of  loads.  Figure  2.4  shows  the  results  of  this  analysis  as  well  as  an  insert 
of  the  beam  itself.  The  results  consist  of  fragmentation  expectation  as  a 
function  of  overpressure  over  a  range  of  physical  size  ranges.  It  is 
interesting  to  note  that  the  results  verify  what  one  would  intuitively  expect: 
that  is,  under  low  load  the  cantilever  beam  will  fail  near  its  support  while 
as  the  pressure  is  increased  the  possibility  of  getting  smaller  and  smaller 
size  particles  becomes  increasingly  more  probable.  Each  particle  size  range 
reaches  a  peak  at  a  particular  load  level  (not  necessarily  predominating). 

At  a  certain  load,  there  is  no  chance  of  getting  any  size  other  than  the 
smallest  size  range  considered.  This  leads  to  the  possibility  of  achieving 
a  pulverization  pressure  if  a  small  enough  size  range  is  considered.  At 
any  specific  load  level  (e.  g.  20  psi)  some  physical  size  intervals  may 
dominate  (e.  g.  8-10  ft,  6-8  ft,  and  0  to  2  ft  )  while  others  may  have  been 
phased  out  or  have  not  been  developed  yet.  (e.  g.  2-4  ft  and  4-6  ft)  This 
is  wholly  dependent  on  the  statistical  nature  of  the  material  (i.  e.  the 
Weibull  parameters  in  this  case  M  =  7,  v ^  =  50,  and  o'  -  1  500). 
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2.  2  Fragment  Trajectory  Analysis 

Trajectory  analysis  for  structural  fragments  starts  with  the 
equation-of-motion  of  a  particle  acted  upon  by  drag  forces.  This  is  a 
second-order  nonlinear  differential  equation,  which  can  be  solved  in  a 
series  of  steps  as  the  fragment  is  followed  through  its  horizontal  translation. 
The  flight  time  of  the  particle  is  determined  by  its  initial  height,  depending 
on  its  location  in  the  structure  under  consideration.  Families  of  trajectories 
can  be  found  for  a  range  of  particle  sizes,  each  set  corresponding  to  some 
combination  of  weapon  yield  and  ground  range.  An  indication  of  debris 
distribution  can  then  be  made  by  combining  all  the  sets  of  trajectories 
corresponding  to  a  specific  weapon  yield. 

The  initial  conditions  required  for  the  use  of  this  transport  model 
are  characterized  by  three  results  of  the  fragmentation  solution,  namely: 
time  of  failure,  size  of  fragment,  and  initial  fragment  velocity.  Solutions 
are  oblaixxed  fur  particle  sizes  consistent  with  those  examined  in  the  frag¬ 
mentation  model  and  with  the  conservative  assumption  that  time  of  failure 
and  failure  velocity  are  zero.  We  now  discuss  the  transport  model  in 
some  detail. 

2.  2.  1  Trajectory  of  a  Particle 

Consider  the  motion  of  a  particle  through  air  such  that  the  drag 
force  acting  on  the  particle  is  proportional  to  the  square  of  the  relative 
'  elocity  between  the  particle  and  the  medium.  It  is  assumed  that  the 
vertical  and  horizontal  motion  of  the  debris  particle  are  uncoupled.  This 
is  true  if  the  center  of  pressure  of  the  particle  coincides  with  its  centroid 
feu*  all  orientations  so  that  no  rotation  occurs,  and  further  that  the  hori¬ 
zontal  component  of  relative  velocity  between  the  air  and  particle  be  signif¬ 
icantly  greater  than  the  vertical  component. 

The  horizontal  equation  of  motion  is  then 

x  -  -----  (x  -  up  '2-6) 

who  re 

x  -  horizontal  coordinate 

y  =  vertical  coordinate,  positive  downward 
(  *  )  =  differentiation  with  respect  to  time 
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c*  =  aerodynamic  coefficient, 
/°  =  mass  density  of  air 
u  =  air  particle  velocity 
g  =  gravitational  constant 

r  ,  >  . 

k  _  r  1;  u  -  x 

1+1;  u  <  x 


(projected  area  x  drag  coefficient) 


mass 


k  =(+  1;  V  £  o 

l;  y  >  o 

This  is  a  Riccati  type  nonlinear  differential  equation  and  can  be  linearized 
by  a  simple  transformation  of  coordinates. 


•  /)  s 

x  =  ft  - 

^  s 


Then 


**07  -0(1' 


Making  the  substitution  into  Eq.  (2-6) 


p 1  -  us'i)2  =  fu2  -  2u  p  i  +  /3  2  (i 2 

^s  ^'vS/  2  '  |  S.  Z  V  s 

The  value  of  ft  can  be  determined  to  make  the  term  vanish. 


Thus 


and  then 


0  =  - 


k  x  /• 


.  2 


.  2 


(2-7) 


(2-8) 


(2-9) 


(2-10) 


*7-0  (!)  =  *0'  u  +2u  !  -0(f) 

ft  —  =  -  0  u2  +  2u  - 

s  s 

ft  ~  s  =  -u^s+  2  ft  u  s 

ft  ^  s  -20us+u^s  =  O 

Treating  u  as  a  constant,  a  closed  form  solution  to  Eq.  (2-11)  can  be  ob¬ 


(2-11) 


tained  in  the  form 


ll 

0  ‘ 


s  =  ( C  .  +  C  t )  e 


(2-12) 
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from  which  follows 


u 


s  =  H  (C]  +c2t)e  ?  * 
where  and  are  constants  of  integration. 
F rom  the  substitution 


+  C2  e 


u  t 


(2-13) 


x  —  (3  - 

s 


It  can  be  shown  that 


x  =  u  + 


C3(3 

T-^t 


(2-14) 


where  the  two  constants  of  integration  C^,  have  been  combined  to  form 
the  single  constant  through  the  relation 


C 


3 


The  horizontal  displacement  follows  by  an  integration  of  Eq.  (2-14)  and  it 
is  given  by 

x  =  c4  +  ut  +  /3  logfi  ( 1  +  C3  t)  (2-15) 

It  will  be  noted  that  a  constant  wind  velocity  u  is  assumed  in  obtaining 
Eq.  (2-14)  and  (2-15).  This  is  not  true  for  the  induced  motion  of  the  particle 
as  a  result  of  blast  winds.  The  wind  velocity,  u,  is  a  function  of  overpres¬ 
sure,  which  is  dependent  on  yield,  ground  range,  and  time.  Hence,  trajec¬ 
tories  must  be  determined  by  computing  the  ho  rizontal  motion  for  a  series 
of  short  time  intervals,  A  t,  over  each  of  which  a  constant  wind  velocity  can 
be  assumed.  Towards  applying  Eq.  (2-14)  and  (2-15)  throughout  the  trajec¬ 
tory  let  us  assume  that  at  the  end  of  the  mth  interval  the  velocity  and  hori¬ 
zontal  displacement  are  x  and  x  .  From  Eq.  (2-14),  setting  t  -  0  one  has 
r  m  m 

C,  =  -i-  (x  -  u  ) 

3  m  8  m  m 
m 

where  u  is  the  (assumed)  constant  wind  velocity  and  {3  is  associated  with 
m  m 

the  (assumed)  constant  local  air  density.  The  remaining  constant  of  integra¬ 
tion  C  is  simply  the  horizontal  displacement  at  the  beginning  of  the  time 
4m 

interval. 
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The  displacement  at  the  end  of  the  time  interval  x  .  ,  is  obtained 
from  Eq,  (2-15)  and  is 

x  ,  ,  =  x  +  u  At+/3  log 
m+1  m  m  m  e 

A  more  directly  useful  form  of  Eq.  (2-16)  is  given  by 

4%=%  A  *  +  loge  (l  +  A  t)  (2-17) 

where  Ax  is  the  increase  in  horizontal  displacement  at  the  time  increment 
m 

At. 


In  Eq.  (2-17),  it  is  assumed  that  u  ,  A  and  x  are  determined  at 
'  m  ~  m  rn 

the  beginning  of  this  interval.  In  sequencing  the  successive  intervals  it  is 
convenient  to  employ  Eq.  (2-14)  to  establish  the  velocity  at  the  end  of  the  time 
interval  At.  Thus,  we  have 


(2-18) 


Consideration  must  also  be  given  to  the  increasing  lag  of  the  particle  behind 
the  shock  front  during  flight,  as  shown  in  Fig.  2.5. 

A  description  of  the  physical  situation  is  as  follows.  Consider  the 
instant  at  which  the  building  is  fragmented.  The  blast  parameters  are  now 
inspected.  Denote  the  weapon  yield  as  W,  ground  range  as  Rq.  Under  these 
conditions  we  can  evaluate  the  shock  wave  velocity,  u,  the  air  mass  density 
f  ,  the  positive  phase  duration  t^,  and  the  overpressure  p  at  the  instant 

the  initial  shock  wave  is  R  distance  from  ground  zero;  or  in  other  words, 

o 

just  as  the  shock  wave  is  exerting  its  influence  on  the  particle.  Now  we 
can  allow  the  action  to  continue  until  some  small  duration  of  time  has  lapsed, 
At.  The  shock  front  hasmoved  a  distance  u  A  t.  From  Eq.  (2-15)  we  can 
find  the  particle  movement,  Xj.  Furthermore,  the  particle  has  fallen 
behind  the  shock  front  because  it  has  not  been  able  to  accelerate  fast  enough 
to  keep  up  with  it.  At  this  point,  we  can  evaluate  a  new  u, p  ,  t,  and  p. 
However,  this  must  be  done  for  the  blast  wave  which  is  now  acting  on  the 
particle.  This  requires  taking  into  account  an  additional  decay  in  overpres¬ 
sure.  After  these  new  "constants"  are  evaluated,  we  can  again  solve  Eq. 
(2-15)  for  some  i.me  increment,  A  t,  and  obtain  the  new  distance  the  particle 
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G  OF  PARTICLE  BEHIND  SHOCK  FRONT 


travels,  x^.  So  far  the  particle  has  traveled  a  total  distance  of  +  x^,  if 
we  do  not  change  At,  the  number  of  times  we  can  increment  is 


V 


2  h 
_ c 

g 


At 


Where 


2  h 


is  the  fall  time  of  a  particle  from  height  h  ,  using  the 


g  *  o 

considerations  outlined  for  each  increment  or  solution  of  Eq.  (2-15),  the 

total  distance  the  particle  travels  is: 


D  = 


i  =  1 


x. 

i 


If  we  index  the  size  of  the  particles  to  show  that  they  are  of  different  sizes, 
we  have  the  total  distance  each  particle  travels  by  repeating  the  solution  of 
Eq.  (2-11)  under  these  conditions,  but  for  different  sizes  of  particles.  We 
obtain: 

n 


v 


n 


x 


2. 

i 


If  we  select  k  sizes: 


n 


It  is  then  possible  to  incorporate  all  this  information  into  graphical  form 
valid  for  the  W  ,  R  in  question,  as  depicted  in  Fig.  2.  6. 

o  o 
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COMPU  I  ER  FLOW  CHART  FOR  PARTICLE  TRAJECTORY  CALCULATIONS 


2.  2.  2  Terminal  Velocities  of  Debris  Particles 


It  is  possible  to  obtain  the  velocities  of  the  different  size  debris 
p?.rticles  in  an  analogous  manner  utilizing  Eq.  (2-11).  Appendix  D  contains 
trajectory  curves  for  debris  transport  that  consider  four  parameters;  dis¬ 
placement,  velocity,  particle  size,  and  starting  height,  for  various  overpres¬ 
sures  due  to  1  MT.  One  outstanding  fact  that  presents  itself  is  the  extremely 
high  terminal  velocities  that  all  particles  have  (e.  g.  ,  in  excess  of  60  ft/sec 
for  a  5  in.  particlt  at  8  psi,  1  MT).  This  bombardment  of  debris  is  bound 
to  have  serious  tertiary  effects  on  both  structural  and  human  elements  in  an 
area  of  high  debris  concentration. 
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CHAPTER  THREE 


DYNAMIC  FLEXURE  EXPERIMENTS  WITH  BRITTLE  BEAMS 


3-  PURPOSE 


In  the  previous  chapter,  a  u  uque  fragmentation  model  was  developed. 
During  the  course  oi  formulating  this  theory  a  fundamental  assumption  was 
made.  This  was  that  the  maximum  dynamic  stresses  introduced  into  the 
various  unit  segments  of  the  loaded  beam  were  independent  of  the  fracture 
characteristics  of  the  segments.  The  experiments  described  below  and 
the  further  analytical  studies  in  Appendix  C  tend  to  verify  the  validity  of 
the  above  assumption. 

Another  purpose  in  conducting  the  following  experiments  was  to  get 
a  better  qualitative  feel  for  brittle  elements  under  dynamic  load.  To  facil¬ 
itate  this,  extensive  high  speed  photography  was  employed  to  record  a  variety 
of  support  and  load  conditions. 

Finally,  an  attempt  was  made  to  investigate  the  fragment  distribution 
of  a  cantilevered  hydrostone  beam  under  a  uniform  dynamic  load. 


3.  *  Static  Versus  Dynamic  Flexure  in  Brittle  Materials 

The  classical  flexure  experiment  with  brittle  beams  is  aimed  at 
establishing  static  bending  strength  of  a  material  (i.  e.  ,  its  modulus  of 
rupture).  Such  tests  are  usually  performed  using  three  or  four  point 
bending  on  a  simply  supported  beam  specimen.  In  the  case  of  three  point 
bending,  for  example,  maximum  flexural  stress  occurs  at  the  center  of  the 
span  for  a  uniform  beam  (i.  e.  ,  under  the  point  of  load),  and  fracture  will 
initiate  there  too  as  the  ultimate  load  is  reached.  This  simple  static  test 


results  in  two  pieces  (fragments),  and  is  consistent  specimen  to  specimen 
with  a  scatter  in  the  magnitude  of  the  failure  stress  as  influenced  by  severity 
of  the  critical  flaws  at  the  center  of  the  span.  In  such  a  case  as  this,  simple 
theory  is  adequate  to  predict  the  location  of  the  failure  section,  but  the  actual 
material  strength  can  only  be  predicted  after  consideration  of  the  statistical 
nature  of  the  material  has  been  made.  Weibull,  for  example,  has  developed 
one  such  statistical  theory  which  can  be  used  to  predict  static  material 


strength  tor  a  variety  of  stress  states  if  certain  "material  parameters’ 


known. 
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When  the  case  of  dynamic  flexure  is  considered  the  problem  becomes 

a  great  deal  more  complicated.  In  the  dynamic  case,  the  elastic  response 

of  the  beam  is  a  function  of  its  dynamic  excitation  (load-time  history);  an 

infinite  number  being  possible.  Thus,  one  begins  by  not  having  a  very  good 

idea  of  what  the  actual  dynamic  bending  moment  diagram  might  be.  Since 

this  quantity  is  time  dependent  and  continually  influenced  by  inertial  forces, 

one  is  confined  to  predicting  stresses  for  relatively  simple  cases  unless 

he  is  prepared  to  undertake  the  numerical  analysis  done  in  Appendix  C. 

When  the  problem  of  predicting  dynamic  stresses  is  combined  with  the 

dynamic  effects  of  statistical  flaw  theory,  prediction  of  the  dynamic  rupture 

characteristics  in  brittle  beams  becomes  very  complicated  indeed.  A 

number  of  investigators  *  have  been  concerned  with  the  dynamic  response 

of  reinforced  concrete  beams,  but  few  have  been  concerned  with  the  dynamic 

rupture  behavior  of  brittle  beams  per  se.  English,  *  *  in  1951,  reported 

a  theoretical  study  of  the  rupture  of  brittle  beams  under  impulsive  loading. 

That  study  was  intended  to  provide  a  pattern  of  the  kind  of  fragmentation 

that  would  result  from  dynamic  flexure,  and  was  based  primarily  on  the  ' 

12  A 

Griffith  Flaw  Theory,  and  the  Weibull  Statistical  Theory  of  Rupture. 

3.  2  Material  Properties 

The  brittle  beams  used  in  this  study  were  cast  of  Hydrostone,  a 
white  cementitious  material  which  is  similar  to  plaster  of  paris.  The 
characteristics  and  some  mechanical  properties  of  the  material  are  given 
below  in  Table  3.  1.  In  this  table,  the  static  and  dynamic  values  of  tensile 
strength  were  measured  using  different  specimen  types  and  methods  of 
test.  Static  tensile  strength  was  measured  by  the  use  of  a  dog-bone  type 
specimen  in  a  direct  pull  test.  Dynamic  tensile  strength  was  measured 
by  the  use  of  a  stress  wave  technique  described  in  detail  in  Reference  13.  In 
this  latter  method  a  traveling  longitudinal  stress  wave  on  a  long  slender 
free -free  bar  caused  tensile  fracture  upon  reflection  of  the  first  compres¬ 
sive  pulse  from  a  free  end.  By  the  use  of  suitably  placed  strain  gages  and 
dynamic  recording,  the  ultimate  fracture  strain  was  monitored  and  converted 
to  ultimate  stress  through  a  modulus  of  elasticity. 

Trade  name  -  U.  S.  Gypsum  Company. 
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Table  3.  1 

CHARACTERISTICS  OF  HYDROSTONE + 


Specific  Gravity 

1. 86  , 

Unit  Weight 

116  lb/ ft 

Static  Tensile  Strength** 

1020  psi 

Static  Modulus  of  Elasticity 

2.  73  x  10^  psi 

+  *f 

Dynamic  Tensile  Strength 

1730  psi 

Dynamic  Modulus  of  Elasticity 

2.  85  x  10^  psi 

111 

Static  Flexural  Strength 

1440  psi 

Static  Compressive  Strength 

6310  psi 

All  material  was  mixed  at  maximum  consistency  with 
0.4  percent  sodate  retarder  added  to  provide  workability. 

+j- 

1/1  inch  square  cross  section 

+++ 

1/2  inch  square  beams  in  pure  bending 

•f  •f  ■f  -r 

3  inch  diameter  cylinders  in  direct  compression. 

The  dynamic  modulus  of  elasticity  used  was  calculated  from  observed 

5 

values  of  the  velocity  of  propagation  of  the  longitudinal  wave  (1.  28  x  10  in. 
per  sec)  and  the  mass  density  per  unit  volume  of  the  material.  The  apparant 
disparity  between  the  static  and  dynamic  values  of  tensile  strength  reported 
in  Table  3.  1  can  be  attributed  to  the  types  of  tests  used.  The  direct  pull 
test  used  for  the  static  tensile  strength  measurements  suffered  from  induced 
bending  caused  by  eccentricity  of  loading.  Hence,  the  measured  static 
tensile  strengths  are  lower  than  they  wo  id  be  if  the  parasitic  bending  effects 
were  eliminated.  The  static  and  dynamic  modulus  of  elasticity  values  agree 
within  6  percent. 

3.  3  Dynamic  Flexure  Experiments 

The  apparatus  used  in  the  dynamic  flexure  experiments  of  this  study 
is  pictured  in  Fig.  3.  1.  The  device  was  essentially  a  steel  frame  with 
hardened  steel  pin  roller  supports  which  prevent  the  development  of  hori¬ 
zontal  reaction  forces.  The  supports  were  placed  5  in.  center  to  center  of 
supports.  The  beams  used  were  about  6  in.  in  length,  and  1/4  in.  square. 
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With  this  arrangement  an  effective  span  to  depth  ratio  of  20:1  was  obtained. 
Experiments  were  performed  by  placing  a  beam  in  the  flexural  device, 
as  shown  in  Fig.  3.  1  and  then  a  charge  of  lead  azide  was  placed  at  the 
center  of  the  span  in  a  small  gelatin  capsule.  The  charge  was  ignited  by 
use  of  a  crimped  nichrome  "hot  wire"  operated  from  a  battery. 

For  some  of  the  experiments  other  end  support  conditions  were 
investigated.  For  these  cases  other  fixtures  were  constructed.  For  the 
case  of  the  "fixed-fixed"  beam,  a  device  providing  clamped  ends  was  used. 

In  the  case  of  the  cantilever  beam,  a  clamp  fixture  at  one  end  of  the  beam 
was  used. 

The  simplest  of  photographic  techniques  was  used  during  the  experi- 
ments.  A  Fastex  high  speed  movie  camera  (approximately  7500  frames  per 
sec)  was  aligned,  sighted,  and  focused  on  the  span  of  the  beam  between  the 
steel  support  frame.  Proper  lighting  behind  the  specimen  through  a  frosted 
glass  plate  provided  a  clear  image  for  photographing.  A  Goose  synchronizing 
unit  was  used  to  coordinate  the  camera  with  the  explosive  impulse.  Initially, 
the  camera  was  started  from  a  switch  and  after  1/2  sec  of  running  time, 
the  explosive  charge  was  detonated  by  the  Goose.  The  entire  experiment 
requires  only  3/4  of  a  second.  The  first  1/2  sec  period  was  used  to  allow 
the  camera  to  get  up  to  speed  and  the  remaining  1/4  sec  was  used  to  detonate 
the  explosive  and  record  the  dynamic  event  on  film. 

Figures  3.  2  through  3.  6  each  show  three  successive  frames  from 
typical  experiments.  In  Fig.  3.  2,  for  example,  the  clear  span  between 
supports  is  plainly  indicated.  The  flag  bearing  the  symbols  "3  SS"  on  the 
left  lower  corner  of  the  frame  indicates  that  a  3  mg  charge  of  lead  azide 
was  used  at  the  center  of  the  span  and  that  the  beam  was  simply  supported. 

The  dark  i  lass  in  the  upper  center  of  the  frame  is  a  cone  shaped  "smoke 
stack"  which  was  used  to  remove  the  smoke  which  would  otherwise  obscure 
the  photographs.  The  vertical  bar  in  the  background  on  the  right  side  of 
tli e  frame  is  a  portion  of  a  ring  stand  which  was  used  to  support  the  "smoke 
stack"  and  the  ignition  wire.  The  single  frames  are  about  150  microsec 
apart. 

Trade  names  -  The  Wollensak  Corporation. 
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The  fracture  which  resulted  in  Fig.  3.  2  was  representative  of  a 
typical  static  break.  One  complete  fracture  occurred,  at  the  center  of 
the  span  (under  the  point  of  loading)  and  only  two  fragments  were  produced. 
This  behavior  was  found  to  be  reproducible  and  typical  of  all  simply  sup¬ 
ported  dynamic  experiments  where  very  small  explosive  charges  were  used. 

The  behavior  which  was  found  to  be  typical  of  the  large  load  levels 
with  simply  supported  beams  is  pictured  in  Fig.  3.  3.  In  this  case  an 
explosive  impulse  of  30  mg  of  lead  azide  was  used.  The  photograph  reveals 
the  extent  of  multiple  fracturing  which  occurred  at  this  high  dynamic  load 
level.  From  a  strictly  qualitative  argument,  the  symmetry  of  the  fractured 
pieces  seems  to  indicated  relatively  simultaneous  formation  and  stochastic 
independence  of  the  failures.  If  this  were  not  true,  an  antisymmetrical  frag¬ 
ment  pattern  would  result. 

Experiments  with  fixed-fixed  beams  showed  results  similar  to  those 
of  the  simply  supported  cases.  At  low  (Fig.  3.4)  load  levels,  three  breaks 
(four  fragments)  were  observed.  This  is  what  again  would  be  expected  in 
the  static  case  of  a  fixed-fixed  beam  since  under  center  point  loading,  this 
beam  will  develop  equal  maximum  bending  moments  at  the  center  of  the 
span  (positive  moment)  and  at  the  end  supports  (negative  moment).  For 
fixed-fixed  beams  having  high  load  levels  multiple  symmetrical  breaks 
nearly  identical  to  those  observed  for  the  simply  supported  beam  were 
found  to  occur.  (Fig.  3.  5) 

One  experiment  was  performed  using  a  cantilevered  beam  and  a 
high  load  level  (30  mg).  In  this  experiment,  one  end  of  the  beam  was 
charged  and  the  charge  was  placed  1  in.  from  the  free  end  (i.  e.  4  in.  from 
the  clamped  support).  Figure  3.  6  illustrates  the  dynamic  behavior  of  the 
brittle  cantilever  for  the  loading  imposed.  In  the  photograph  the  free  end 
of  the  beam  is  on  the  left  side  of  the  frame.  The  figure  dramatically  shows 
that  localized  fracture  and  multiple  fragments  are  Produced  near  the  free 
end  of  the  beam  where  the  dynamic  load  was  imposed. 

Since  the  result  shown  in  Fig.  3.  6  is  contrary  to  static  analysis,  a 
dynamic  analysis  was  made  with  a  lumped  mass  model  as  outlined  in  Ap¬ 
pendix  C.  The  results  of  this  analysis  confirm  the  behavior  observed. 
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Fig.  3,2  SIMPLY  SUPPORTED  BEAM  UNDER  CENTER  POINT 
DYNAMIC  LOADING  -  3  mg  LOAD  LEVEL 


3.  3  SIMPLY  SUPPORTED  BEAM  UNDER  CENTER  POINT 
DYNAMIC  LOADING  -  30  mg  LOAD  LEVEL 


Fig.  3.  4  FIXED  ENDED  BEAM  UNDER  CENTER  POINT  DYNAMIC 
LOADING  -  5  mg  LOAD  LEVEL 
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Fig.  3.  5  FIXED  ENDED  BEAM  UNDER  CENTER  POINT  DYNAMIC 
LOADING  -  30  mg  LOAD  LEVEL 


Fig.  3.  6  CANTILEVER  BEAM  UNDER  END  POINT  DYNAMIC 
LOADING  -  30  mg  LOAD  LEVEL 


It  is  interesting  to  note  that  for  the  load  imposed  the  maximum  allowable 
bending  moment  was  exceeded  in  0.1  microsec.  Since  the  velocity  of  propaga¬ 
tion  was  observed  to  be  1.  28  x  10"*  in.  /sec  and  the  masses  were  lumped 
1/2  in.  apart,  the  wave  takes  at  least  6  microsec  to  reach  the  nearest  adjacent 
masses.  These  masses,  however,  have  long  since  exceeded  their  allowable 
maximum  in  bending.  Thus  it  is  now  shown  in  a  quantitative  as  well  as 
qualitative  manner  that  fracture  may  develop  simultaneously  and  independently 
at  more  than  one  station  along  a  brittle  beam.  This  effect  was  shown  for 
a  concentrated  loading.  In  the  more  practical  case  of  uniform  dynamic 
loading,  this  phenomena  would  be  even  more  pronounced;  not  -mg  con¬ 
strained  to  a  localized  region. 

3. 4  Shock  Tube  Experiments 

Shock  tube  experiments  were  carried  out  to  determine  the  piece  size 
distribution  versus  pressure  relationship  for  fixed  ended  hydrostone  beams 
under  uniform  dynamic  pressure  loadings.  The  experimental  scheme  is 
illustrated  in  Fig.  3.  7.  Beams  1/4  in.  square  and  6  in.  long  were  mounted 
in  a  fixture  which  attaches  to  the  open  end  of  the  shock  tube.  The  fixture 
clamps  both  ends  of  the  beams,  providing  a  clear  span  of  5  in.  on  a  "fixed- 
fixed”  beam.  A  burst  diaphragm  containing  an  explosive  charge  at  its 
center  is  installed  some  distance  upstream  at  the  static  load  section  of  the 
tube.  The  explosive  charge  is  then  hooked  into  the  igniter  circuit.  Next, 
a  static  pressure,  P^,  is  placed  in  the  static  load  section  from  an  ordinary 
nitrogen  gas  bottle.  When  the  tap  switch  in  the  igniter  circuit  is  closed, 
the  explosive  charge  bursts  the  diaphragm  and  the  resulting  shock  wave 
travels  down  the  tube  and  breaks  the  beams  dynamically.  An  open-end  shock 
tube  was  used  to  avoid  the  effects  of  wave  reflection  ordinarily  encountered 
in  closed  ended  shock  tubes. 

Dynamic  pressures  were  measured  by  the  use  of  piezoelectric  pres¬ 
sure  gages  placed  at  the  center  and  end  of  span.  The  output  of  these  gages 
was  displayed  on  a  suitable  oscilloscope,  equipped  with  a  scope-back  Polar¬ 
oid  camera  which  made  the  dynamic  recordings.  Results  of  the  dynamic 
pressure  measurements  showed  very  little  difference  between  the  pressures 
at  midspan  and  at  the  fixed  support.  Based  on  this  result  the  pressures 
measured  in  each  experiment  were  averaged  and  taken  as  the  uniform 
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dynamic  pressure  over  the  clear  span.  Static  loading  pressures  ranging 
from  zero  to  92  psi  where  used  to  charge  the  tube  and  the  resulting  average 
downstream  pressures  on  the  beams  varied  from  13.  5  to  31.  8  psi.  The 
lowest  pressures  obtained  resulted  from  the  use  of  an  explosive  charge  only, 
with  no  pressure  in  the  static  load  section  of  the  tube.  Since  only  about 
8  psi  would  be  required  to  cause  static  failure  in  the  beams,  the  minimum 
dynamic  pressure  loading  greatly  exceeds  the  pressure  required  to  cause 
failure. 

Results  in  the  range  13.  5  to  31.  8  psi  tend  to  confirm  our  intuitive 
and  analytical  observation  that  piecesize  decreases  with  increasing  load. 
This  is  shown  in  Fig.  3.  8  which  relates  the  mean  fracture  length  for  a 
ten  unit  beam  to  increasing  load. 
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Fig.  3.8  MEAN  FRACTURE  LENGTH  VERSUS  PRESSURE 


CHAPTER  FOUR 

TYPICAL  STRUCTURAL  APPLICATION  OF  DEBRIS  MODELS 


4.  INTRODUCTION  . 

This  chapter  deals  with  the  application  of  the  methods  developed  in 
the  previous  chapter  to  predict  how  debris  resulting  from  a  large  multistory 
structure  forms  and  is  transported  offsHe. 

4.  1  Model  Structure 

The  model  structure  to  be  analyzed  is  a  40  story  curtain  wall.  The 
walls  are  faced  with  brick  which  is  6  in.  thick.  The  total  wall  height  is 
340  ft  and  the  panel  height  is  8-1/2  ft.  The  support  structure  of  the  wall 
will  be  considered  infinitely  ductile;  thus,  only  the  outside  brick  material 
in  the  panels  will  be  considered  as  possible  transportable  debris. 

Figure  4.  1  illustrates  the  wall  described  above  except  that  window 
openings  will  not  be  considered  in  the  fragmentation  analysis  of  an  individ¬ 
ual  panel. 

4.  2  Loading 

It  is  anticipated  that  within  the  fireball  region  of  an  atomic  blast, 
total  destruction  will  occur.  For  situations  such  as  this,  debris  distri¬ 
bution  may  be  made  by  the  methods  developed  in  the  previous  report.  * 

This  report  is  concerned  with  those  regions  which  are  characterized  by 
only  partial  destruction.  This  partial  distribution  may  include  failure 
of  exterior  and  interior  walls,  but  excludes  the  total  collapse  and  transition 
on  an  entire  structure  (i.  e.  floor  panels,  support  structure,  etc.).  The 
loading  that  will  be  investigated  will  therefore  be  from  about  8  psi  to  an 
overpressure  region  just  short  of  total  destruction  (18  psi). 

Figure  4.  2  indicates  the  blast  wave  approach  to  the  frontal  face 
of  the  structure.  It  is  seen  that  within  the  overpressure  region  of  interest 
the  structures  will  be  well  within  the  range  of  the  mach  stem.  It  is  therefore 
felt  that  the  pressure  will  be  uniformly  distributed  across  the  frontal  face 
of  the  building.  Since  each  bay  on  each  story  will  contain  one  window#  the 
building  should  undergo  diffraction  type  loading  exclusively,  at  least  for 
the  walls. 
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The  frontal  wall  may  be  expected  to  receive  an  overpressure  in 
excess  of  twice  the  sideon  pressure.  This  phenomenon  occurs  due  to  the 
fact  that  as  the  front  of  the  blast  wave  strikes  the  face  of  the  structure 
reflection  takes  place.  As  the  blast  wave  moves  forward,  the  front  face 
reduces  to  the  sideon  pressure.  The  rear*  or  leeward  face,  of  the 
structure  will  experience  a  pressure  in  the  opposite  direction  to  the  pres* 
sure  on  the  frontal  face.  This  pressure  will  have  a  magnitude  equal  to  the 
sideon  overpressure.  What  then  is  the  effective  loading  on  the  structure? 
This  is  depicted  in  Fig.  4.  3  which  illustrates  the  front  wall  of  the  structure 
fragmenting  due  to  a  load  of  twice  the  sideon  pressure.  The  rear  wall  is 
fragmented  by  a  load  equal  to  the  sideon  overpressure.  Although  the  shock 
wave  caused  fragmentation  in  opposite  directions  on  the  front  and  rear  walls, 
the  pressure  acting  on  the  leeward  wall*  from  the  rear,  comes  off  when  the 
shock  front  is  oast  the  front  wall.  The  only  force  acting  on  the  fragments 
wi*'.  be  due  tc,  le  blast  winds  associated  with  the  passing  shock  front. 

Debris  resulting  from  the  fragmentation  of  the  side  walls  will  not  be  in¬ 
cluded  at  this  time. 

4.  3  Panel  Fragmentation 

Figure  4  4  is  a  beam  analogy  for  the  fragmentation  of  a  masonry 
panel  (i.  e.  ,  the  panel  is  bounded  by  the  enclosing  floors  and  is  between 
two  columns).  Failure  of  a  masonry  wall  panel  has  been  shown  to  be 
characterized  by  the  "arching  effect"  along  most  of  its  length  (except 
near  the  column  supports).  By  expanding  the  fragmentation  distribution 
of  a  rigid  beam  along  this  "arching  length"  most  of  (at  some  load  levels 
all)  the  panel's  masonry  material  will  be  accounted  for.  In  those  cases 
where  the  "arching  effect"  only  acts  for  some  distance  which  is  short 
of  the  panel's  total  length,  the  remaining  material  could  be  characterized 
by  another  beam,  cantilevered,  on  each  of  the  ends  adjacent  to  the  columns, 
and  this  beam's  fragmentation  distribution  could  be  expanded  between  the 
enclosing  floors.  The  panel  to  be  considered  in  the  following  application 
will  have  its  fragmentation  analysis  based  on  the  entire  panel  material  acting 
as  a  r;gid  beam  of  length  equal  to  the  panel  height  (i.  e.  8  1/2  ft)  and  cross 
section  equivalent  to  the  thickness  squared  (i.e.  ,  f>  in,  by  6  in.  ).  The  re¬ 
sults  of  this  analysis  of  the  masonry  rigid  beam  is  exhibited  in  Fig.  4.  5  and 
is  analogous  to  the  results  depicted  in  Fig.  2.4  for  the  cantilevered  hydro- 
stone  beam. 
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Fig.  4.  4  BEAM  ANALOGY  FOR  A  MASONRY  PANEL  Fv\ILURE 
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Fig.  4.  5  FRAGMENTATION  CURVE  FOR  MASONRY  PANEL 


4.  4  Orientation  of  Structure  with  Respect  to  the  Shock  Front 

The  initial  distribution  studies  will  not  consider  any  orientations 
other  than  the  frontal  wall  being  incident  to  the  shock  wave.  It  is  under¬ 
stood  that  orientation  will  play  a  large  role  in  determining  debris  distri¬ 
butions.  from  both  a  loading  and  transport  standpoint.  This  report  will 
not,  however,  include  this  effect. 

4.  5  Fragmentation  Results 

Figure  4.  5  is  the  fragmentation  curve  for  the  masonry  panel 

analyzed  as  a  rigid  masonry  beam  and  discussed  in  Section  4.  3  above. 

The  number  of  trials,  n,  used  in  the  fragmentation  analysis  was  again  15. 

The  curve  indicates  that  fragmentation  occurs  at  about  8  psi  and  total 

pulverization  (i.  e.  only  the  smallest  size  interval,  0-2.  32  in.  ,  is  present) 

is  achieved  at  about  20  psi.  Weibull  parameters  for  the  masonry  beam 

wei  e  determined  by  experience  consistent  with  the  strength  of  the  material. 

(i.  e.  ,  m  =  6,  o'  =  448.0,  and  ^  =  381.0) 
o  u 

4.  6  Transport  of  Fragment  Sizes 

The  curves  in  Fig.  4.6  describe  the  mode,  panel's  fragmentation 
and  transport  characteristics  at  a  10  psi  overpressure  range.  It  is  thus 
possible  to  find  the  expectation  of  obtaining  a  particular  size  particle  a 
given  distance  from  its  preblast  position. 

4.  7  Debris  Distribution  Profiles 

The  information  supplied  in  Fig.  4.6  is  sufficient,  along  with  some 
knowledge  of  the  wall  dimensions,  to  develop  a  debris  profile.  The  pro¬ 
cedure  employed  is  to  uniformly  distribute  the  material  from  each  story 
height  over  its  transported  distance.  For  example,  at  10  psi  the  particles 
of  size  range  2.  32  -  4.  64  in.  equivalent  spherical  radius  at  the  40th  story 
travel  345  to  515  ft  horizontally  from  their  original  position  before  striking 
ground.  We  thus  distribute  32  percent  (i.  e.  the  percentage  of  material  be¬ 
tween  2.  32  and  4.  64  in.  radius  at  10  psi)  of  the  wa.T  material  for  one  panel 
over  an  area  of  length  170  ft  and  unit  width.  The  height  of  this  distributed 
material  may  be  determined  by  dividing  the  percentage  of  fragments  of  size 
2.  32  -  4.  64  in.  radius  by  the  unit  width  area  over  which  the  material  dis- 
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The  above  procedure  for  determining  debris  height  implies  a  unit  volume 
of  material  at  the  40th story  level (i.e.,  the  panel  at  this  level  is  1  ft  long 
by  1ft  thick).  Thus  the  debris  height  is  "normalized"for  a  unit  panel  volume, 
however,  the  actual  thickness  of  the  panel,  utilized  later  in  determining  ac - 
tual  debris  heights,  is  implied  (i.  e.,  6  in. )  implicitly  by  the  development  of 
the  Weibull  parameters  used  in  the  fragmentation  analysis.  If  we  repeat 
the  procedure  for  all  size  ranges  at  a  given  story  and  then  at  each  story 
of  the  wall,  then  by  superimposing  all  these  individual  debris  depths  upon 
one  another  we  obtain  a  normalized  debris  profile.  The  profiles  are 
normalized  in  the  sense  that  we  must  multiply  the  debris  depth  by  the 
volume  of  material  in  the  panel  at  one  story  level,  (i.  e.  ,  8  1/2  ft  high 
by  1/2  ft  thick  by  1  ft  wide  minus  a  correction  for  window  openings).  The 
super  position  procedure  is  illustrated  in  Fig.  4.  7. 

Appendix  E  contains  profiles  for  a  range  of  wall  heights  and  in 
some  cases  fragmentation  criteria  were  varied  over  a  constant  wall  height 
to  illustrate  the  sensitivity  of  the  debris  profile  to  fragment  particle  size. 
The  three  cases  of  fragmentation  criteria  were: 

a)  100  percent  small  particles  0-2.  32  in.  in  radius. 

b)  The  panel  fragmentation  behavior  at  10  psi  as  calculated 
by  the  methods  developed  in  Chapter  II  and  illustrated 

in  Fig.  4.  5  for  the  masonry  beam. 

c)  100  percent  large  particles  9.-28  -  11.  60  in.  in  radius. 

The  critical  characteristics  of  these  curves  are  specified  in  Table  4.  1. 

5|C 

The  parameters  X  and  y  are  the  same  as  illustrated  in  Fig.  1.  3.  Again, 
in  order  to  relate  the  y  to  actual  debris  height  it  must  be  multiplyed  by 
the  unit  width  volume  of  panel  material  at  one  story  level,  (i.  e.  8  1/2  ft 
by  1/2  ft  corrected  for  any  wall  aperatures). 
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Table  4.  1 


SENSITIVITY  OF  DEBRIS  CONTOURS  TO  BUILDING  HEIGHT 

AND  PARTICLE  SIZE 


All  Small 

Predicted 

All  Large 

Building 

Height 

X 

❖ 

y 

X 

❖ 

y 

X 

* 

y 

40 

514 

0.  003 

200 

0.  036 

200 

0.  085 

35 

496 

0.  0027 

180 

0.  035 

180 

0.  080 

30 

460 

0.  0023 

165 

0.  032 

130 

0.  075 

25 

425 

0.  019 

155 

0.  029 

-  - 

-  -  - 

20 

400 

0.  015 

135 

0.  026 

130 

0.  060 

15 

350 

0.  012 

125 

0.  023 

125 

0.  055 

10 

300 

0.  009 

75 

0.  018 

90 

0.  045 

5 

175 

0.006 

50 

0.  015 

40 

0.  040 

CHAPTER  FIVE 


CONCLUSIONS  AND  RECOMMENDATIONS 

5.  SUMMARY 

The  work  in  the  first  four  chapters  of  this  report  completes  this 
phase  of  the  debris  study.  It  is  significant  to  note  that  it  is  now  possible 
to  predict  the  deposition  of  offsite  debris  with  a  one-to-one  correspondence 
between  the  original  structural  material  and  its  final  rubble  form  for  brit¬ 
tle  wall  elements.  This  debris  model  is  based  upon  a  fragmentation  theory 
that  is  derived  from  a  statistical  approach  to  the  strength  of  brittle  materials. 
The  validity  of  the  theory  has  been  verified  by  the  dynamic  experiments 
of  Chapter  Three.  The  fragmentation  model  was  outlined  for  a  variety  of 
end  conditions  and  is  a  significant  computational  improvement  over  the 
method  previously  developed.  * 

Debris  profiles  of  hypothetical  masonry  wall  panels  were  generated 
utilizing  a  simple  distribution  procedure.  The  results  of  this  study  indicate 
the  validity  of  the  general  procedure  for  obtaining  the  debris  profile  for  a 
single  building. 

5.  1  Conclusions 

1.  Debris  may  be  divided  into  two  categories,  offsite  and  onsite. 
Onsite  debris  prevents  no  impedance  to  logistic  activity  other  than  localized 
shelter  rescue  operations. 

2.  In  order  to  characterize  offsite  debris  profiles,  the  piecesize 
distribution  of  the  fragments  must  be  known.  This  is  in  order  to  character¬ 
ize  the  particle's  trajectory. 

3.  Low  level  dwellings  must  be  included  as  potential  debris  producing 
structures.  Even  elements  at  10  ft  above  ground  may  be  transported  signif¬ 
icant  distances  (e.  g.  into  the  street)  at  moderate  overpressure  levels 

(10  psi). 

4.  In  general,  smaller  particles  from  higher  initial  heights  will  be 
transported  greater  distances  than  larger  particles  at  lower  heights.  As 
overpressure  increases  the  size  of  the  particle  fragmenting  will  be  smaller. 
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5.  Particles  fragmenting  and  being  transported  offsite  from  their 
original  position  will  have  an  extremely  high  terminal  velocity  (i.e.,  at  least 
50  ft  per  sec).  For  this  reason,  offsite  debris  must  be  considered  an  extremely 
dangerous  secondary  effect  of  blast.  This  effect  must  be  evaluated  in  light 
of  its  damage  producing  capabilities  to  structures  as  well  as  a  casualty  pro¬ 
ducing  mechanism  to  unsheltered  populence. 

6.  The  assumption  that  the  maximum  dynamic  stresses  introduced 
into  the  various  unit  segments  are  independent  of  the  fracture  character¬ 
istics  has  been  shown  both  experimentally  (Chapter  Three)  and  analytically 
(Appendix  C).  This  establishes  the  theoretical  development  of  the  frag¬ 
mentation  model. 

7.  The  statistical  strength  characteristics  of  typical  unit  segments 
have  been  expressed  by  their  cumulative  distribution  function,  F(p).  This 
function  gives  the  probability  of  fracturing  the  unit  at  a  load  magnitude  equal 
to  or  less  than  p  and  for  computational  purposes  has  been  expressed  in 
Weibull  form.  It  is  to  be  emphasized  that  this  form  was  chosen  for  con¬ 
venience.  At  present,  WeibuT  parameters  are  not  known  for  many  common 
construction  materials  and  must  be  determined  as  outlined  in  Appendix  B. 

8.  Debris  profiles  from  multistory  walls  may  be  generated  as  out¬ 
lined  in  Chapter  Four.  As  one  may  see  from  Appendix  E  the  profiles  are 
typically  unimodal  in  shape  and  of  course  are  very  sensitive  to  piecesize 
distribution.  For  walls  it  5,  10,  15  stories  in  height  profiles  exist  for 
each  of  three  piecesize  distributions,  all  large,  predicted,  and  all  small. 

It  is  interesting  to  note  that  the  large  particles  generally  lend  most  of  the 
shape  to  the  predicted  size  profile.  That  is,  up  to  the  peak  the  large  size 
particles  predominate.  This,  of  course,  is  in  part  due  to  the  method  of 
developing  these  profiles.  Large  particles  are  distributed  over  much  less 
area  than  small  ones  and  hence  tend  to  produce  much  greater  depths. 

9.  The  profiles  in  Appendix  E.  i:ave  their  ordinate  normalized  with 
respect  to  volume  and  packing.  (Packing  corresponds  to  a  void  ratio  of 

l.  0  in  Appendix  E. )  If  the  ordinate  were  multiplied  by  the  unit  width  volume 
of  material  at  one  story  height  and  then  again  by  a  suitable  void  ratio,  then 
the  profile  wouid  express  the  true  depth  of  the  building's  transported  debris 
at  ail  points  along  its  transported  distance. 
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5.  2  Recommendations 


After  reading  this  report  one  can  see  that  the  balance  covers  only 
one  area  as  outlined  in  Chapter  One  (i.  e.,  brittle  wall  elements).  In  addition, 
the  analytical  methods  of  fragmentation  and  subsequent  transport  did  not 
consider  either  the  effects  of  orientation  of  the  structure  to  the  blast  or 
shielding  of  one  building  by  another.  It  is  therefore  recommended  that  any 
follow  on  study  should: 

1.  Continue  to  develop  analytical  methods  for  additional 
structural  materials  and  elements  as  well  as  complete 
structures.  These  analytical  models  should  be  suitably 
verified  by  appropriate  experimental  investigation. 

2.  The  ir  dividual  building  debris  contours  must  be  com¬ 
bined  to  give  a  cumulative  debris  contour  for  an  entire 
subarea  of  contiguous  structures. 

3.  Effects  of  blast  orientation  aud  structural  shielding 
should  be  accounted  for  in  specifying  the  load  on  the 
structure. 

4.  Finally,  all  methods  developed  should  be  combined  to 
yield  a  single  computational  model,  as  outlined  in 
Chapter  One. 
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APPENDIX  A 

COMPUTATIONS  OF  STATIC  STRESS  DISTRIBUTIONS 
FOR  DIFFERENT  END  CONDITIONS 
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Fig.  A.  1  CANTILEVER  BEAM  UNIFORMLY  DISTRIBUTED  LOAD 
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Fig.  A.  2  SIMPLY  SUPPORTED  BEAM-UN1FORMLY  DISTRIBUTED  LOAD 


APPENDIX  B 

DETERMINATION  OF  STATISTICAL  STRUCTURAL  PARAMETERS 


The  statistical  strength  characteristics  of  a  typical  unit  segment  are 
described  by  its  cumulative  distribution  function,  F(P).  This  function  gives 
the  probability  of  fracturing  the  unit  at  a  load  magnitude  equal  to  or  less  than 
the  applied  loading.  The  distribution  function  may  be  determined  by  physi¬ 
cally  testing  many  unit  segments  or  by  appealing  to  one  of  the  "weakest  link" 
statistical  fracture  theories.  Of  these,  Weibull's  is  the  most  papular  theory 
now  in  use. 


One  method  for  determining  the  values  of  Weibull's  statistical  param¬ 
eters  m,  <t  and  nr^  is  the  graphical  method.  The  method  can  be  applied  for 
any  type  of  test,  however,  for  illustrative  purposes  a  rectangular  beam  under 
pure  bending  will  be  considered  here. 


Let  N  be  the  number  of  nominally  identical  specimens  whose  fracture 
strengths  have  been  determined.  The  fracture  strengths  are  ordered  in  in¬ 
creasing  magnitude  such  that 
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Corresponding  to  the  fracture  strength  <rn  of  rank  n,  the  cumulative  probabil¬ 
ity  of  failure  is  given  by 
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For  a  rectangular  beam  under  pure  bending,  Weibull's  expression  for 
the  cumulative  probability  of  failure  is,  for  cr>(r  , 
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here  (r  is  the  maximum  fiber  stress.  This  equation  may  be  rewritten  as 
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It  can  be  seen  that  a  plot  of  this  distribution  function  will  be  linear  in  a  sys^ 
tern  of  rectangular  coordinates  in  which 
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is  tic  ordinate  and  log  (o  -  )  is  the  abscissa. 
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Also,  m+1  will  he  the  slope  of  the  distribution  in  these  coordinates. 

Thus  the  determination  of  the  Wei  bull  parameters  proceeds  as  follows. 
A  guess  is  made  for  the  value  of  <r  .  A  plot  of 

iog  log  |-y-~F — |+log  versus  log  {<rn  -  rj 

is  made  and  fitted  with  a  curve.  If  the  curve  is  concave  up,  the  value  of  cr^ 
should  be  decreased  (however,  a  negative  <ru  is  physically  unreasonable)  and 
if  the  curve  is  concave  down,  the  value  of  cr^  should  be  increased.  This  pro  - 
cess  is  pursued  until  the  resulting  curve  is  (or  may  be  approximated  by)  a 
straight  line.  This  value  of  cr^  is  the  correct  one.  The  value  of  m  is  deter¬ 
mined  as  the  slope  minus  one.  Taking  V  to  be  the  total  volume  of  the  beam 
in  the  same  units  of  length  as  appearing  in  the  stresses  a-  and  cr^,  the  value 
of  or  may  be  determined  by  inserting  the  coordinates  of  a  point  on  the  line 
into  Eq.  (4). 
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APPENDIX  C 


A  DYNAMIC  ANALYSIS  OF  A  CANTILEVER 


UNDER  AN  IMPULSIVE  LOAD 


For  the  beam  shown  below  a  concentrated  load  P  is  applied  at  an 
arbitrary  distance  (a)  from  the  free  end. 

P..X, 

M 


b  =  i  -a 

a 

1  l 

Fig.  C.  1  CANTILEVER  UNDER  ARBITRARY  CONCENTRATED  LOAD 


From  elementary  beam  theory  the  deflection  x  distance  from  the  free  end  is: 

y  =  -44t-?-  (3b-i  +  x)  (C-l) 


We  now  let  f. .  represent  the  deflection  at  point  x.  due  to  a  unit  load  (Pel) 

J 

at  point  x..  Then  using  Eq.  (C-l)  and  letting  (a)  represent  x.,  x  replace 
J  J 

x.  and  P  =  1  2 

f..  =  ...urn.1—-  (3  ( i  -  x  ?  -  l  +  x  ),  i>j  (C-2) 
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Formula  (C-2)  generates  the  lower  triangular  portion  of  the  flexibility 
matrix  [>]  •  Let  the  cantilever  in  Fig.  C-l  be  approximated  by  n  equally 

lumped  masses  as  shown  below  and  loaded  at 
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Fie.  C.2  LUMPED  MASS  APPROXIMATION 
mass  number  p.  Then  each  lumped  mass  is 
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where/5  is  the  beam  mass  pi  r  unit  length. 

1 1 T  RESEARCH  INSTITUTE 


-6'»- 


The  equilibrium  equations  which  represent  Fig.  C.  2  are 


Thus 
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Equation  (C-3)  represents  a  system  of  second  order  differential  equations 
with  the  load  P  at  mass  p  a  function  of  time.  This  problem  is  an  initial 
value  one  with  the  following  conditions: 

Y  =  0  at  t  =  0 

and  Y  =  0  at  t  =  0 

The  initial  value  problem  was  solved  using  a  Runge-Kutta  numerical 
integration  scheme.  At  any  time,  t,the  moment  at  point  i,  M.,  is 
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Figure  C.  3  indicates  the  load-time  history  corresponding  to  a  30  mg  charge 
which  was  used  for  (t).  Figure  C.4  shows  the  nondimensional  bending 
moment  diagram  M/M  versus  x  for  various  t. 

°  0  u  a2 

i  o’  B  k  “ 

where  M  = - , - 

o  6 

and  o-  is  the  dynamic  tensile  stress. 
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Distance  from  Free  Fud  of  Beam,  Indies 
Fig.  C.  4  DYNAMIC  BEAM  BENDING  MOMENT 
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APPENDIX  D 


FRAGMENTATION  AND  TRANSPORT  CURVES  FOR  W  =  1  MT 

Trajectory  calculations  relating  distance  traveled  to  impact  velocity 
for  various  initial  heights  were  carried  out,  according  to  the  method 
described  in  Section  2.  2.  The  procedure  utilized  equivalent  spherical 
fragment  sizes  corresponding  to  the  results  of  the  masonry  fragmentation 
study  of  Chapter  Four.  Spherical  fragments  were  selected  for  this  analysis 
since  drag  effects  are  identical  for  all  orientations  of  the  sphere;  whereas 
little  is  known  of  drag  effects  of  irregular  material  at  random  orientation*, 
or  during  rotation.  The  effect  of  this  assumption  is  to  average  out  the 
effects  of  large  and  small  projected  areas  from  the  irregular  objects. 

A  value  of-  1.  0  was  used  for  the  drag  coefficient,  C  ..  This  is  an 
approximate  value  in  the  velocity  range  the  particles  travel.  The  selected 
value  of  the  material  density  was  135  lb  per  ft  ,  which  is  in  the  range  of 
most  masonry  building  materials. 

Results  of  the  calculations  performed  on  the  IBM  7094/1401  digital 
computer  are  plotted  in  the  following  figures. 
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Fii*.  D.  I  TP  '  NS  PORT  CURVES  AT  11,400  ft  FROM  1  MT  SURFACE  BURST 


i'OH  I  crr,V!-.S  Ai’  10,400  ft  FROM  1  MT  S U R FACE  BURST 


Fie.  D.  3  TRANSPORT  CURVES  AT  9,  900  ft  FROM  1  MT  SURFACE  BURST 
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l  1  TRANSPORT  CURVES  AT  7,  400  ft  FROM  1  MT  SURFACE  BURST 


APPENDIX  E 

DEBRIS  PROFILES  FOR  HYPOTHETICAL  WALL;  W  =  1  MT 

Debris  profiles  were  developed  for  the  example  wall  as  outlined  in 
Chapter  F our.  In  addition  to  varying  the  height  of  the  wall,  the  fragment 
piecesize  distribution  was  also  varied  in  some  cases. 
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I  Fig.  E.  2  DEBRIS  PROFILE  FOR  A  WALL  HEIGHT  OF  ^  STORIES 

i  EXPOSED  TO  10  psi  AT  I  MT 


I 


-89- 


Distance  from  Wail,  ft  (x  10  ) 


Fig.  E.  3  DEBRIS  PROFILE  FOR  A  WALL  HEIGHT  OF  30  STORIES 
EXPOSED  TO  10  psi  AT  1  MT 
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Fig.  E.  5  DEBRIS  PROFILE  FOR  A  WALL  HEIGHT  OF  25  STORIES 
EXPOSED  TO  10  psi  AT  1  MT 
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Distance  from  Wail,  ft  (x  10  ) 


Fig.  E.  6  DEBRIS  PROFILE  FOR  A  WALL  HEIGHT  OF  20  STORIES 
EXPOSED  TO  10  psi  AT  1  MT  (All  Small  Particles) 
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DEBPIS  PROFILE  FOR  A  WALL  HEIGHT  OF  15  STORIES 
EXPOSED  TO  10  psi  AT  1  MT  (ALL  Large  Particles) 
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Fig.  E.  1  1  DEBRIS  PROFILE  FOR  A  WALL  HEIGHT  OF  10  STORIES 
EXPOSED  TO  10  psi  AT  l  MT  (All  Small  Particles) 
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Normalized  Debris  Height,  ft  (x  10  ) 
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Fig.  E.  1  1  DEBPIS  PROFILE  FOR  A  WALL  HEIGHT  OF  10  STORIES 
EXPOSED  TO  10  psi  AT  1  MT  (All  Large  Particles) 
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